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. . . Abstract 
00 

' Considering analyticity, unitarity and an influence of coupled channels, ex- 

perimental data on the isovector P-wave of vrvr scattering was analyzed to study 
yo-like mesons below 1900 MeV. The analysis indicates evidently that in the 
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^ I energy range 1200-1800 MeV, there are three p-like mesons: /3(1250), /9(1450) 

and /?(1600), unlike the PDG tables. The obtained P-wave Trvr-scattering length 
(N ! (a} = 33.9 ± 2.02 [lO-^m-^*]) most matches to the one calculated in the local 

Nambu-Jona-Lasinio model. 

Oh 1 Introduction 

I— I, The investigation of vector mesons is actual up to now due to their role in forming 
■ the electromagnetic structure of particles and because, e.g., in the p-family, only the 
>■ p(770) meson can be deemed to be well understood [1]. The other p-like mesons must be 
either still confirmed in various experiments and analyses or their parameters essentially 
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corrected. For example, the p(1250) meson which was discussed actively some time 
. ago [2, 3] was confirmed relatively recently in the amplitude analysis of the LASS 
^ ' Collaboration [4] and in combined analysis of several processes [5]. However this state 
[ — is referred to only slightly in the PDG issue [1] (the relevant observations are listed 
O ■ under the p(1450)). 

'i> ■ On the other hand, the vrvr interaction plays a central role in physics of strongly 

^ . interacting particles and, therefore, it has always been an object of continuous inves- 

^ ' tigation. Let us note only some recent works devoted also to the theoretical study of 
■ " " ' the isovector P-wave of tttt scattering. First, there are the analyses of available experi- 
mental data on the vrvr scattering utilizing the Roy equations [6, 7, 8] and the forward 
dispersion relations [9, 10], in which, e.g., low-energy parameters of the vrvr scattering 
were obtained. Second, there are the works in which the low-energy parameters are 
calculated in chiral theories with the linear realization of chiral symmetry [11, 12]. 

We use our model-independent method [13]. It is based on the first principles 
(analyticity and unitarity) directly applied to analysis of experimental data. Our aim 
is to study the p-like mesons below 1900 MeV and to obtain the vrvr-scattering length. 
Unfortunately this method, using essentially an uniformizing variable, is applicable 
only to the 2-channel case and under some conditions to the 3-channel one. Here the 
thresholds of the vrvr and uir channels are allowed for explicitly in the uniformizing 
variable (in the threshold region of the second channel, one has observed a deviation 
from elasticity of the P-wave vrvr scattering). Influence of other coupled channels is 
supposed to be taken into account through the background. In order to investigate a 
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coupling of resonances with these other channels, we also apply multichannel Breit- 
Wigner forms to generate the resonance poles. 

The paper is organized as follows. In Section II, we outline the method of the 
uniformizing variable in applying it to studying the 2-channel nn scattering and present 
results of the analysis of the available data [14]-[16] on the isovcctor P-wave of tttt 
scattering. Section III is devoted to analysis of the same data using the Breit-Wigner 
forms. Finally, in Section VI, we summarize and discuss obtained results. 

2 Analysis of P-wave nn scattering in the uniformiz- 
ing variable method 

Let the TTTr-scattering j^-matrix be determined on the 4-sheeted Riemann surface with 

the right-hand branch-points at 4m^+ and (m^^ + 771^^0)"^ and also with the left-hand one 
at s = 0. It is supposed that influence of other branch points can be taken into account 
through the background. The Riemann-surface sheets are numbered according to the 
signs of analytic continuations of the channel momenta 



as follows: signs(lmA;i, Im/i;2) = -|— 1-, — h, and H — correspond to sheets I, II, III 

and IV, respectively. 

The (S'-matrix is supposed to be S — SresSbg where Sres represents resonances and 
Sbg, the background. In general, an explicit allowance for the (m^j-l-mTro)^ branch point 
would permit us to describe transitions between the vrvr and wn initial and flnal states 
with the help of the only one function d{ki,k2) (the Jost matrix determinant) using 
the Le Couteur-Newton relations [17]. Unfortunately, the data on process tttt —>■ ujtt 
are absent. 

In Ref. [13] it was shown how one can obtain the multichannel resonance repre- 
sentations by poles and zeros on the Riemann surface with the help of the formulas, 
expressing analytic continuations of the matrix elements, describing the coupled pro- 
cesses, to unphysical sheets in terms of those on sheet I. It is convenient to start from 
resonance zeros on sheet I. Then in the 2-channel tttt scattering, we have three types of 
resonances: (a) described by a pair of complex conjugate zeros in the (S-matrix element 
on sheet I and by a pair of conjugate shifted zeros on sheet IV; (b) described by a pair 
of conjugate zeros on sheet III and by a pair of conjugate shifted zeros on sheet IV; (c) 
which correspond to a pair of conjugate zeros on sheet I, the one on sheet III and two 
pairs of conjugate zeros on sheet IV. Due to unitarity the poles on sheet II, III and IV 
are situated in the same energy points as the corresponding zeros on sheet I, IV and 
III, respectively. Note that the size of shift of zeros on sheet IV relative to the ones on 
sheets I and III is determined by the strength of coupling of the channels (here tttt and 
lot:). The cluster kind is related to the nature of resonance. 




and 
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With the help of the uniformizing variable^ 



(nr. + /;/-^(i)/2 \Js — liir^ , + y .s — + iiIt^o)'^ 



^7 {{m^ + m^o)/2f -ml+ 



the considered 4-shceted Riemann surface is mapped onto the ?;-plane, divided into two 
parts by a unit circle centered at the origin. Sheets 1 (11), 111 (IV) are mapped onto the 
exterior (interior) of the unit disk on the upper and lower f -half-plane, respectively. The 
physical region extends from the point i on the imaginary axis (tttt threshold) along the 
unit circle clockwise in the 1st quadrant to the point 1 on the real axis (a;7r° threshold) 
and then along the real axis to the point b — {m^ -\- 171^^0 + 2m^+)/ {m^^ + m^o — 2m7r+) 
into which s = 00 is mapped on the v-plane. The intervals (—00, —b], [—b~^, b~^], [b, 00) 
on the real axis are the images of the corresponding edges of the left-hand cut of the 
TTTT-scattering amplitude. The (a) resonance is represented in S'(7r7r tttt) by two pairs 
of poles on the images of sheets 11 and 111, symmetric to each other with respect to the 
imaginary axis, and by zeros, symmetric to these poles with respect to the unit circle. 
Note that a symmetry of the zeros and poles with respect to the imaginary axis appears 
due to the real analyticity of the ^-matrix, and the symmetry of the poles and zeros 
with respect to the unit circle ensures a realization of the known experimental fact that 
the TTTT interaction is practically elastic up to a vicinity of the ujn^ threshold. 

The resonance part of S'-matrix Sres becomes a one-valued function on the I'-plane 
and, in the tttt channel, it is expressed through the d(i')-function as follows^: 

Sres — ^7/ N ^ (2) 

d{v) 

where d{v) represents the contribution of resonances, described by one of three types 
of the pole clusters in the 2-channel case, i.e., 

M 

d = v-^'l[il-v:v){l + Vr,v) (3) 

n=l 

with M the number of pairs of the conjugate zeros. 
The background part S^g is taken in the form: 



Sbg = exp 



y (ai + a2 e{s - si) + as e{s - S2) 



(4) 



^The analogous uniformizing variable has been used, e.g., in Ref. [18] in studying the forward elastic 
pp scattering amplitude and in Ref. [19] in the combined analysis of data on processes tttt — > tttt, KK 
in the channel with I'^J^^' = 0+0++. 

^ Other authors also have used the parameterizations with the Jost functions at analyzing the S- 
wave WW scattering in the one-channel approach [20] and in the two-channel one [21]. In latter work, 
the uniformizing variable ^2 has been used and the TTTr-thrcshold branch-point has been neglected, 
therefore, their approach cannot be employed near by the ww threshold. 
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where = aj + ibi, si is the threshold of Air channel noticeable in the p-like meson 
decays, S2 is the threshold of p27r channel. Due to allowing for the left-hand branch- 
point at s = in the uniformizing variable (1), Oi = 61 = 0. Furthermore, 62 = is an 

experimental fact. 

With formulas (2)-(4), we have analyzed data [14]-[16] for the inelasticity parameter 
(77) and phase shift of the TTTr-scattering amplitude (5) (S'(7r7r nn) = r] exp{2id)), 
introducing three (p(770), p(1250 - 1580) and p(1550 - 1780)), four (the indicated 
ones plus p(1860 — 1910)) and five (the indicated four plus p(1450)) resonances. In 
Rcf.[16], results of two analyses arc cited: one uses the s-channel helicity amplitude 
when extracting the Trvr-scattcring amplitude on the vr-cxchangc pole; in the other, the 
t-channel one is used instead. Therefore, we have taken both analyses as independent. 
There are given the data for the phase shift of amplitude below the KK threshold. 
Comparing these data with the ones of Refs.[14, 15], one can see that the points of 
the former lie systematically by l°-5° higher than the ones of the latter, except for 
two points of Ref. [14] at 710 and 730 MeV, which lie by about 2° higher than the 
corresponding points of Ref. [16] and which are omitted in the subsequent analyses. 
Since we do not know the energy dependence of the remarked deviations of points, we 
have supposed a constant systematic error that must be determined in the combined 
analysis of data. We have obtained a satisfactory description with x^/NDF and with 
the indicated systematic error equal respectively to 293.518/(186 — 15) = 1.716 and 
-1.885° for the case of three resonances, 280.043/(186 - 19) = 1.677 and -1.897° for 
four resonances, 269.574/(186 — 23) = 1.654 and —1.876° for five resonances. When 
calculating for the inelasticity parameter, three points of data [15] at 990, 1506 
and 1825 MeV have been omitted in all three cases as giving the anomalously big 
contribution to x^- When calculating foi' the phase shift, three points of data [16] 
have been omitted in all three cases: the one at 790 MeV from the s-channel analysis, 
and two at 790 and 850 MeV from the i-channel one. 

On figures 1, we demonstrate results from our fitting to data [14]-[16]. The short- 
dashed curves are for the three-resonance description, the long-dashed ones for the 
four-resonance and the solid for the five-resonance cases. 

From a variety of possible resonance representations by pole-clusters, the analyses 
prefer he following one to be the most relevant: the p(770) is described by the cluster 
of type (a.) and the others by type (b). 

Let us show the obtained pole clusters on the lower -^/s-half-plane in Table 1 (it is 
clear there are also complex-conjugate poles on the upper half-plane). The background 
parameters arc: 02 = 0.0093, 03 = 0.0618 and 63 = —0.0135 for the thrcc-rcsonance, 

02 = 0.00166, 03 = 0.0433 and 63 = -0.00442 for the four-resonance and 02 = 0.0248, 

03 = 0.0841 and 63 = 0.0019 for the five-resonance descriptions. 

Though the description can be considered practically as the same in all three cases, 
careful consideration of the obtained parameters and energy dependence of the fitted 
quantities suggests that the resonance p(1900) strongly desired and that the p(1450) 
should be also present improving slightly the description. 

Masses and widths of the obtained p-states can be calculated from the pole positions 
on sheet II for resonances of type (a) and on sheet IV for resonances of type (b). If the 
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7T + 71 ^ 71 + 7T" 
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Figure 1: The phase shift of amphtude (the upper figure) and module of matrix ele- 
ment (the lower figure) of the P-wave TTTr-scattering. The results of 3-resonance (short- 
dashed), 4- resonance (long-dashed), and 5- resonance (solid) descriptions are shown. 
The data from Refs. [14]-[16]. 



resonance part of the amplitude reads as 

rpres "\/^ -^e; 

^les -S-iy/s Ftot ' 

we obtain for the masses and total widths, respectively, the following values (in the 
MeV units): 



for 


p(770), 


769.3 


and 


146.6; 


for 


p(1250). 


1256.2 


and 


261.4; 


for 


p(1470). 


1468.8 


and 


200.8; 


for 


p(1600). 


1595.1 


and 


145.8; 


for 


p(1900). 


1895.4 


and 


175.2. 
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Table 1: Pole clusters for the p-like resonances. 



Three resonances 


Sheet 


II 


III 


IV 


p(770) 


MeV 


767.3 - i73.3 


782 - i65.6 




p(1250) 


MeV 




1219.9 - /■152 


1219 - /■116.2 


p(1600) 


MeV 




1585 - il30.5 


1578 - i72.2 


Four resonances 


Sheet 


II 


III 


IV 


p(770) 


MeV 


766.5 - i73.2 


783.1 - i66.2 




p(1250) 


MeV 




1251.4 -il52.1 


1249 - il44.3 


p(1600) 


MeV 




1585.2 -il41. 8 


1579.6 - i73.6 


p(1900) 


MeV 




1871.5 -i97.2 


1894 - i95.3 


Five resonances 


Sheet 


II 


III 


IV 


p(770) 


ys;, MeV 


765.8 - i73.3 


778.2 - i68.9 




p(1250) 


MeV 




1250 - il31.4 


1249.4 -il30.7 


p(1470) 


MeV 




1469.2 - i89.3 


1465.4 - ilOO.4 


p(1600) 


MeV 




1634.8 - il45.9 


1593.4 - i72.9 


p(1900) 


MeV 




1883 - il06.5 


1893.4 - i87.6 



3 The Breit— Wigner analysis of P-wave nn scatter- 
ing 

In various works [1], it was shown that the p-like resonances obtained in the previous 
section have also other considerable decay channels in addition to those considered 
explicitly above. It was observed that the p(1450) and/or a possible p(1250) can also 
decay to the rjp^ (< 4%), (pTi (< 1%), and An (seen) channels, where the fraction Fj/F 
[1] is given in parenthesis. The p(1700) resonance has the large branching ratio to the 
4:71 (large), p27r (dominant) and rjp^ (seen) channels. 

To include explicitly influence of some selected channels and to obtain information 
about couplings with these channels on the basis of analysis of the TTTr-scattcring data, 
we have used 5-channel Breit-Wigner forms in constructing the Jost matrix determi- 
nant d{ki, ■ ■ ■ , k^) . To generate the resonance poles and zeros in the ^'-matrix, the 
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Le Couteur-Newton relation was utilized: 



d{-ki, ■ ■ ■, fcs 
d{ki,---,k5) 



(5) 



In eq.(5) ki, k2, k^, /C4 and kr, are the Tin-, 7r+7r~27r°-, 2(7r+7r^)-, ri27i- and cjvr'^-channcl 
momenta, respectively. The ci-function is taken as 0? = dresdbg with dres, describing 
resonance contributions, and d^g, the background. 

The Breit-Wigner form for the resonance part of the d- function is assumed as 



dres{s) = n 



,f2 



(6) 



where Prj — ki{s)/ki{M^), f^j/Mr is the partial width of resonance with mass M^, Rrj 
is a Blatt-Weisskopf barrier factor [22] conditioned by the resonance spins. For the 
vector particle this factor have the form: 



R 



M2 - 4m2 rrjf 



(7) 



with radius r^j = 0.7035 fm which in our analysis are equal for all resonances in all 
channels. Furthermore, we have assumed that the widths of resonance decays to the 
7r"'"7r~27r° and 2(7r"'"7r~) channels are related by relation: fr2 — frs/V^- 
The background part d^g is 



d, 



bg = exp 



Oil + a2 



Sl 



(8) 



where ctj = a, + ibi, Si is the threshold of p27r channel (it is clear that 02 and 62 take 
into account also influence of other channels opened at higher energies than the p27r 
threshold); bi is taken to be zero. 

Using formulas (5)-(8) we have performed the analysis, just as in the previous 
section, with three, four and five resonances. We have obtained the same reason- 
able description in all three cases: the total xV^DF = 316.206/(186 - 17) = 1.871, 
314.688/(186 - 22) = 1.919, and 303.101/(186 - 27) = 1.906 for the case of three, 
four, and five resonances, respectively. On the figure 2, we demonstrate results of our 
fitting to data [14]- [16] for the case of five resonances. The systematic error of data 
[16], discussed in the previous section, is equal to —1.987° in this case. When calculat- 
ing for the inelasticity parameter, three points of data [15] at 990, 1030 and 1825 
MeV were omitted as giving the anomalously large contribution to x^. When calcu- 
lating for the phase shift, the same three points of data have been omitted as in 
the model-independent analysis. For the background we find: Oi = —0.00121 ± 0.0018, 
02 = -0.1005 ± 0.011 and 62 = 0.0012 ± 0.006. 

The obtained values of resonance parameters for the last case are given in Table 2. 
Note that in Table 2 we do not show the errors for /42, /43, /44 and /45, because the 
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7T + 71 ^ 7T" + 7T" 
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Figure 2: The phase shift of amphtude and module of matrix element of the P-wave 
TTTT-scattering. The curves show result of fitting to data [14]-[16] using the Breit-Wigner 
form. 

present data are not sufficient to fix reliably these parameters. We may conclude that, 
using the accessible data, the model-independent analysis testifies in favour of existence 
of the p(1900) and, maybe, p(1450), but the Breit-Wigner approach cannot verify this 
result. 

We have also calculated the isovector P-wave length of tttt scattering a\. Its value 
is shown in Table 3 in comparison with the results from various evaluations in the local 
[12] and non-local [11] Nambu-Jona-Lasinio (NJL) model and from the ones with the 
use of Roy's equations [9, 8, 7] (in Ref.[8] one apply also the chiral perturbation theory 
(ChPT) to construct a precise TTTT-scattering amplitude at s^^'^ < 0.8 GeV). 

4 Conclusions 

The reasonable description of all the accessible experimental data on the isovector P- 
wave of TTTT scattering for the inelasticity parameter (j]) and phase shift of amplitude (5) 
[14]- [16] have been obtained up to 1880 MeV based on the first principles (analyticity 
and unitarity) directly applied to analysis of the data. Analysis has been carried out in 
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TaV)le 2: The /-j-like resonance parameters (all in MeV). 



State 


p(770) 


p(1250) 


p(1450) 


p(1600) 


p(1900) 


M 


777.69±0.32 


1249.8±15.6 


1449.9±12.2 


1587.3±4.5 


1897.8±38 


frl 


343.8±0.73 


87.7±7.4 


56.9±5.4 


248.2±5.2 


47.3±12 


frl 


21.6±5.8 


186.3±39.9 


100.1±18.7 


210.2±8.6 


73.7 


fr3 


34.8±8.2 


263.5±56.5 


141.6±26.5 


339.7±12.5 


104.3 


frA 




231.8±111 


141.2±98 


141.8±33 


9 


frb 




231±115 


150±95 


108.6±40.4 


10 


^tot 


Ril54.3 


>175 


>52 


>168 


>10 



Table 3: Comparison of values of the rnr scattering length a\ from various approaches. 



a\[lQ-'m-l] 


References 


Remarks 


33.9 ±2.02 


This paper 


Breit-Wigner analysis 


34 


[12] 


Local NJL model 


37 


[11] 


Non-local NJL model 


37.9 ±0.5 


[8] 


Roy equations using ChPT 


38.4 ±0.8 


[9] 


Roy equations 


39.6 ±2.4 


[7] 


Roy equations 



the model-independent approach using the uniformizing variable (here the satisfactory 
description is obtained: x^/NDF = 1.654) and applying multichannel Breit-Wigner 
forms to generate the resonance poles and zeros in the S'-matrix (x^/NDF = 1.906). The 
aim of analysis (except for obtaining an unified formula for the P-wavc tttt scattering 
amplitude in the whole of investigated energy range) was to study the p-like mesons 
below 1900 MeV and to obtain the P-wave TTTr-scattering length. 

For the p(770), the obtained value for mass is a little smaller in the model-independent 
approach (769.3 MeV) and a little bigger in the Breit-Wigner one (777.69 ±0.32 MeV) 
than the averaged mass (775.5 ± 0.4 MeV) cited in the PDG tables [1]. However, this 
mass also occurs in analysis of some reactions [1]. The obtained value of the total width 
in the first analysis (146.6 MeV) coincides with the averaged PDG one (146.4±1.1 MeV) 
and it is a little bit bigger in the second analysis (ps; 154.3 MeV) than the PDG value, 
however, this is also encountered in other analyses [1]. Note that predicted widths of 
the p(770) decays to the 47r-modes are significantly larger than, e.g., those evaluated in 
the chiral model of some mesons based on the hidden local symmetry added with the 
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anomalous terms [23]. 

The 2nd p-like meson has the mass 1256.2 MeV in the 1st analysis and 1249.8±15.6 
MeV in the 2nd one. This differs significantly from the mass (1459±11 MeV) of the 2nd 
p-like meson cited in the PDG tables [1]. We mentioned already in Introduction that 
the p(1250) meson was discussed actively some time ago [2, 3], and next the evidence 
for it was obtained in some analyses [4, 5]. Note also the talk by Ichiro Yamauchi at 
the HADRON'07 conference [24] in which this mass value of the second p-like meson 
is also confirmed in the re-analysis of some data on the e+e~ — > cutt^ reaction. To the 
point, if this state is interpreted as a first radial excitation of the l+l state, then 
it lies down well on the corresponding linear trajectory with an universal slope on the 
(n, M^) plane [25] (n is the radial quantum number of the qq state), while the meson 
with mass M — 1450 MeV turns out to be considerably higher than this trajectory. 

Further we note that some time ago it was shown that the 1600-MeV region contains 
in fact two p-like mesons. This conclusion was made on the basis of investigation of 
the consistency of the 27r and Ati electromagnetic form factors and the Trvr-scattering 
length [26] and as a result of combined analysis of data on the 2tt and in final states 
in e'^e~ annihilation and photoproduction [27]. We assume this possibility, i.e., that in 
the energy range 1200-1800 MeV there are three p-like mesons and the neoct meson has 
the mass about 1450 MeV. This does not contradict to data; in the model-independent 
analysis, description is even slightly better than without this state. It seems that this 
improvement of description can become more remarkable at an explicit consideration 
of the 7727r-threshold in the uniformizing variable, because the Breit-Wigner analysis 
points to the considerable coupling of the p-like mesons with the r]2Tr channel. How- 
ever for the final conclusion the combined analysis of several processes to which the 
investigated resonances are appreciably coupled have to be performed. 

The fourth p-like meson turns out to have the mass about 1600 MeV rather than 
1720 Me F cited in the PDG tables [1]. Note that in a number of previous analyses of 
some reactions one has also obtained the resonance with mass near 1600 MeV [1] . Note 
also a rather big obtained coupling of these p-like mesons with the An channels. 

As to the p(1900), in this region there are practically no data on the P-wave of 
nn scattering. The model-independent analysis, maybe, somehow testifies in favour of 
existence of this state, whereas the Breit-Wigner one gives the same description with 
and without the p(1900). For more definite conclusion about this state, the P-wave nn 
scattering data above 1880 MeV are needed. Furthermore, the combined analysis of 
coupled processes should be carried out. 

Finally, we note that the P-wave nn scattering length {a\ = 33.9 ±2.02 [10~^m~+]), 
obtained in the Breit-Wigner analysis, matches most to the one in the local Nambu- 
Jona-Lasinio model [12]. 
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